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Abstract 

We present the domain-wall fermion operator which is reflection symmetric in the fifth dimension, 
with the approximate sign function S{H) of the effective 4-dimensional Dirac operator satisfying the 
bound |1 —S'(A)| < 2dz for G ^max]^ where dz is the maximum deviation |1 —-v/x-Rz(2^)|max 

of the Zolotarev optimal rational polynomial Rz{x) of Ijy/x for x G [1, with degrees 

(n — 1, n) for Ng = 2n, and (n, n) for Ng = 2n + 1. 


1 


I. INTRODUCTION 


In lattice QCD with exact chiral symmetry [l-3|, the overlap Dirac operator with bare 
quark mass in general can be written as 


D{mg) =mg + 


(1 — rrriq) 


2r 


1 + 75 


H 

y/JP, 


, r = l/[2mo(l - hmo)], mo G (0,2), (1) 


where H = cih^(l + and c and d are constants. Here = ■y^Dn,, and is the 

standard Wilson-Dirac operator minus the parameter mo G (0,2). 

The eigenmodes of ([I]) consist of complex conjugate pairs, and (for topologically non¬ 
trivial gauge held) real eigenmodes with dehnite chiralities at rUg and 1/r satisfying the 
chirality sum rule jdj, — n_ -|- iV_|_ — N_ = 0, where {N±) denote the number of 
eigenmodes at rUg (1/r) with ± chirality. Empirically, the real eigenmodes always satisfy 
either (n_ = = 0, n_|_ = N_) or (n_|_ = N_ = 0, n_ = A^+). Thus, we have 


det D{mg) = 


{rmgY+ defH^ = {rrrig) "+ defH^, > 0, 

/2 _ ^ > 0 , 


(rmg)"“ det = (rm^) 


where = P±{D^D), and P± = (1 ± 75)/2. It follows that the pseudofermion action 
for any number of havors of overlap fermion can be expressed in terms of n± and 
(Hermitmn and positive-dehnite), thus is amenable to the hybrid Monte Carlo simulation 
(HMC) 


fl. 


However, this approach requires the computation of the change of n± at each 
step of the molecular dynamics in HMC, which is prohibitively expensive for large lattices 
0, B]- Moreover, the discontinuity of the fermion determinant at the topological boundary 
highly suppresses the crossing rate between different topological sectors, thus renders HMC 
failing to sample all topological sectors ergodically. 

These difficulties can be circumvented by using domain-wall fermion (DWF) with hnite 
Ng in the hfth dimension. Then HMC of lattice QCD with DWF on the 5-dimensional 
lattice can sample all topological sectors ergodically and also keep J^h^ chiral symmetry at 
a gooc 
QCD 


precision. This has been demonstrated for 2-flavors QCD (8|, [Oj, and (l-l-l)-flavors 


10l |. For DWF with hnite Ng, it is vital to preserve the chiral symmetry maximally. 


or equivalently, the approximate sign function S{H) of the effective 4-dimensional Dirac 
operator satishes the bound, |1 —S'(A)| < dz for G [AP„, A^^^,], where dz is the maximum 
deviation |1 — y/xRz{x)\uiax of the Zolotarev optimal rational approximation Rz{x) of l/\/x 
for X G [1, A^^^/A^j„], with degrees (n, n) for Ng = 2n + 1, and {n — 1, n) for Ng = 2n. This 
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can be attained by assigning a weight to each layer (along the hfth dimension) of the DWF, 


according to the formula 


111 


OJ., = 


Xr 


-y/l - K'^sn^ k'), s = !,■■■, Ns, 


( 2 ) 


where sn(ns; k') is the Jacobian elliptic function with argument Vg (see Eq. (13) of Ref. Q) 
and modulus k' = — 1/6, b = X^^^^/XN^. Here XN^ are the upper-bound and 

lower-bound for the eigenvalues of H^. It should be emphasized that Xmin and Xmax have to 
be hxed properly for each set of simulations, depending on the parameters [3 = quark 
masses, and lattice size, such that the desired precision of chiral symmetry can be attained 
with the minimal cost of the simulation. 

Nevertheless, ([2]) breaks the (reflection) symmetry in the hfth dimension, which is 
essential for obtaining the exact pseudofermion action for hybrid Monte Carlo simulation 
of one-havor DWF 10|], as well as other applications. In this paper, we obtain the weights 
satisfying the symmetry (ws = uns-s+i, 5 = 1 , • • • , Ng) such that the approximate sign 
function S{H) of the effective 4-dimensionaI Dirac operator satishes the bound, |1 — S'(A)| < 
2dz-i which is twice of that of the optimal DWF without R^ symmetry [n|. In other words, 
if one imposes the R^ symmetry, this is the optimal chiral symmety one can have. 


II. OPTIMAL DOMAIN-WALL FERMION 


In general, the 5-dimensionaI lattice operator of DWF can be written as 



\3-Xij^')\xx'',gg' (^Ps^w R ^')xx'bgs’ {(^gDyj 1 ) 3 . 3 ;/(3) 

where pg = cujg + d, ag = cug — d, and c, d are constants. The indices x and x' denote the sites 
on the 4-dimensionaI space-time lattice, and s and s' the indices in the hfth dimension, while 
the lattice spacing a and the Dirac and color indices have been suppressed. The operator L 
is independent of the gauge held, and it can be written as 


L — P+L^ + P-L_, — (1 ± 75 ) 72 , 


and 


{L+')ss' 


(T—)s's 


Tfldjsf^^gi, s 1 , 

Js— 1 ,S ^5 1 7 5 ^ Ng, 
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where Ng is the number of sites in the hfth dimension, m = rniq, niq is the bare quark 
mass, and r = l/[2mo(l — drrio)]. Including the action of the Pauli-Villars fields (with bare 
mass mpv = 1/^), the partition function of DWF in a gauge background can be integrated 
successively to obtain the fermion determinant of the effective 4-dimensional Dirac operator. 


Z = 



exp = det D{mq), 


where 


D{mq) = {Dc + mq){l + rDc) = d- 

^ 1 1 + 7^3 (H) 


(1 


rm„ 


2r 






S{H) = 


= H = cH^{l + d'y 5 H^)-\ 


(4) 


Hw (h) 


i + nSi^; i+u;gH‘ 

For the optimal DWF without symmetry [Ul, the weights {cu^} are fixed according 
to the formula ([2]), then S{H) is equal to the Zolotarev optimal rational approximation of 
the sign function of H, i.e, S{H) = HRz{H‘^), satisfying the bound, |1 — S'(A)| < dz for 
G [^miny Knax]y wliere dz is the maximum deviation |1 — ■\/xRz{x)\raax of the Zolotarev 
optimal rational approximation Rz{x) of 1/y/x for x G [1, A^„,,,/A^j„], with degrees (n, n) 
for Ng = 2n -|- 1, and (n — 1,?7,) for Ng = 2n. However, it breaks the R^ symmetry, i.e., 
R^V{rn)R^ ^ V{m), where (i? 5 )s,<i' = <5s',Vs+i-s- In fhe following, we construct the optimal 
DWF satisfying the R 5 symmetry. 


III. OPTIMAL DOMAIN-WALL FERMION WITH R 5 SYMMETRY 


First we recall the basic features of the optimal rational approximation of the ratio of 
two positive and continuous functions f{x)/g{x) for x G [!,&]• Let R^^’'^\x) denote an 


irreducible rational polynomial of the form 


PmX"^ + Pm-lX"^ ^ H - h PO 

QnX^ + qn-ix^~^ H-h go 


{n > m, pi.Qi > 0). 


According to de la Vallee-Poussin’s theorem and Chebycheff’s theorem, the necessary and 
sufficient condition for R^'^’'^\x) to be the optimal rational approximation of f{x)/g{x) for 
X G [1,6] is that 5{x) = f{x) — (x) has (m -|- n -|- 2 ) alternate change of sign in the 

interval [ 1 , 6 ], and attains its maxima and minima (all with equal magnitude), say. 


6(a:) = +A, -A, • • • , (-l)"^+"+iA 
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(a) (b) 

FIG. 1: The deviation 6{x) = 1 — of the Zolotarev optimal rational polynomial 

R^^’^\x) for X G [1,^], b = ^maxl^min ~ (6.20/0.05)^ = 15376. In the left figure (a), {m,n) = 
{n — l,n) = (7,8), while in the right hgure (b), {m,n) = (n,n) = (8,8). In both cases, d{x) has 
(m + n + 2) alternate change of sign in the interval [1,6], with (m +1) minima and (n +1) maxima, 
all of the same magnitude, satisfying the necessary and sufficient condition of the optimal rational 
approximation of Ijy/x for x G [1,6], acccording to de la Vallee-Poussin’s theorem and Chebycheff’s 
theorem. 


at consecutive points (xj, z = 1, • • • , m + n + 2) 

1 = Xi < 0:2 < ■ ■ • < Xm+n+2 = b. 

For f{x)/g{x) = Ij^/x, Zolotarev obtained two optimal rational polynomials in 1877, 


R 


(n—l,n) 


x) and R 


{n,n) 


X] 


Q, i 


in terms of the Jacobian elliptic functions. The basic for¬ 


mulas of Zolotarev optimal rational polynomial are collected in the Appendix. (See also 
Refs. 14, 15].) 

For the optimal DWF without R^ symmetry [h], the approximate sign function S{H) 
of the effective 4-dimensional lattice Dirac operator is exactly equal to HR^^~^’^\h‘^) for 
Ns = 2n, and HR^z’'^\h‘^) for Ns = 2n + l. 

In Fig. [H we plot the deviation = 1 — y/xR^P’"'\x) of the Zolotarev optimal 

rational polynomials Rz’^\^) Rz’^'^ for x G [1,6], 6 = ^max/^min = (6.20/0.05)^ = 

15376. In both cases, 6{x) has (m -|- n -|- 2) alternate change of sign in the interval [1,6], 
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with (m + 1) minima and (n + 1) maxima, all of the same magnitude, thus satisfying the 
necessary and sufficient condition for the optimal rational approximation of l/\/x in the 
interval [1, b], acccording to de la Vallee-Poussin’s theorem and Chebycheff’s theorem. Note 
that dz — 1.19447x 10“^ in the left hgure (a), which is about 2.2 times of dz — 5.39351 x 10“® 
in the right figure (b), i.e., ~ 2.2di”’"'\ as pointed out in Ref. 15|. 


A. Ng = 2n (even) 


For the case Ng = 2n (even), the requirement of symmetry implies that Ug = cuvs+i-s, 


and 


(5(A) = 1 - ^(A) = 


2 nr=i(i “ 


> 0 , 


( 6 ) 


n:=i(i-a;.A)2+m=i(i+^sA)2 

unlike the (5^ ^’"’^(x) (shown in Fig. [I](a)) which varies alternatively between +dz and 
—dz- Nevertheless, one can shift 6 z{x) by a constant +dz, simply by changing the overall 


coefficient Dq of flTS]) to Dq, 


p/ _ ^ nr=i('^+ C21-1) 


, 2^^2 ( K '. ,\ ■ 


K sn^ 


2n 


;k') 


(7) 


Then the maximal deviation becomes 2dz, which is twice of that of the Zoloterev optimal 
rational polynomial. Obviously, this is the optimal chiral symmetry one can have for DWF 
with i ?5 symmetry. This can be attained by requiring d(A) to have 2n+l alternative maxima 
and minima for A G [A^m, >^max], 

1 - A 


(5(A) — 2dz, 0, ■ • ■ , 2dz, 0, 2dz, dz — 


( 8 ) 


1 + A’ 

at the consecutive points (A^m = Ai < A 2 < • • • < A 2 n+i = ^max), where dz is the maximum 
deviation |1 — y/xR^^ ^’"^^(3^)|max of the Zolotarev optimal rational polynomial 
for l/\/x, X G [1,5]- This immediately implies that the weights in ([3]) can be hxed by the 
positions of the minima of S{x), 


^Ns+l—s 




(9) 


where sn(M; k') is the Jacobian elliptic function with modulus td = — 1/6, and K' is 

the complete elliptic function of the first kind with modulus k'. Then the approximate sign 
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30x10-® 


25x10-® 



0.05 0.1 0.2 0.5 1 2 4 6.2 

A, 

FIG. 2: The deviation 1 — 5(A) of the optimal DWF with symmetry, for Ng = 16, Xmaxl^min = 
6.20/0.05, where the weights are computed according to Q. 


function S{H) of the effective 4-diniensional Dirac operator (jl]) becomes 


S(H) = ACidh-' 


h = 


H 


( 10 ) 


which is exactly the same as the Zolotarev optimal rational approximation of the sign func 
tion (see Eq. flT8|) ) except replacing Dq by Dq. Note that if we replace Dq by Dq, 

nr=i(i+ ^ 21 -^ 


Do = 


Y[ViY + C 2 i) ’ 


( 11 ) 


then 5{x) is shifted by a constant —dz, thus its maxima become 0, and minima —2dz- 
However, this cannot be realized by DWF with even Ng and R 5 symmetry, due to the 
constraint (]6]), (5(x) > 0. 

In Fig. 121 we plot 1 — S{x) of the optimal DWF with symmetry, for Ng = 16, 
^max/^min = 6.20/0.05, and the weights according to ([9]). We see that 1 — S{x) has 2n + 1 = 
17 alternate maxima and minima in the interval [Xmin, Amax], with 8 minima and 9 maxima, 
and |1 — 5(A)|maa: = 2(7^, where dz is the maximum deviation |1 — y/xR^^’^'^\max of the 
Zolotarev optimal rational polynomial, as shown in Fig. [11(a). 
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B. Ng = 2n + l (odd) 


For the case Ng = 2n + 1, the requirement of symmetry implies that ojg = ojns+i-s-, 
and 


(5(A) = 1 - ^(A) 


2(1 — Ci;„_|_iA) nr=i(^ “ 

(1 —cu^+iA) ns=i(i ~ + (1 +w^+iA) ns=i(i+’ 


which can be positive, zero, or negative. Thus, we have two options to attain the optimal 
chiral symmetry with symmetry, i.e., to shift (5(A) by a constant Rdz, or —dz- 

The hrst option is to require 5(A) to have 2n + 2 alternative maxima and minima in the 
interval [A^^m) A^q^,], 


5(A) = 24, 0, • • • , 24, 0, 4 


1 — a 
1 + (t’ 


at consecutive points {Xmin = Ai < A 2 < • • • < A 2 n +2 = X^ax), where dz is the maximum 
deviation |1 — ’”'^(x)|max of the Zolotarev optimal rational polynomial R^^'"'\x) for 

1/4^ in the interval [1,5]. Thus the weights in (|3]) can be hxed by the positions of the 
minima of 5(a;), 


— ^A^s + 1 —S 



K'^sn^ 


/(2s- l)iF' 

V 'Ns 



s = 1,- 


, {Ng + l)/2. 


( 12 ) 


Note that the smallest weight oJn+i is chosen to be the unpaired one at the center of the 
hfth dimension, such that 5(A) does not change sign for all A G [Xmin, Xmax]- Then the 
approximate sign function S{H) of the effective 4-dimensional Dirac operator (|1]) becomes 


S{H) 


. Utiih^ + C2i) 
^^UtYh^ + C2i-iy 


h 


H 

^min 


(13) 


which is exactly the same as the Zolotarev optimal rational approximation of the sign func¬ 
tion (see Eq. flT^ ) except for the overall coefficent (ig. 



J_ TT (/ + C 2 I- 1 ) 
\// (^ + C 21 ) 




K'^sn^ 




(14) 


The second option is to require 5(A) to have 2n + 2 alternative maxima and minima in 
the interval [Xmin, Xmax], 


5(A) = 0 , - 24 , • • ■ , 0 , - 24 , 4 
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1 — a 
1 + a' 
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(a) 

FIG. 3: The deviation 1 
and Xrnax 

m- 





(b) 

S{X) of the optimal DWF with symmetry, for Ng = 17, Xmin = 0.05, 
6.20. In (a), the weights are computed according to (fT^ . while in (b), according to 


at consecutive points (Xmin = Xi < X 2 < ■ ■ ■ < \ 2 n +2 = ^max)- This immediately implies 
that the weights in ([3]) can be hxed by the positions of the maxima of S{x), 


UJg UJ — g 


1 

^min 



K'^sn^ 


f (Ng + 3-2s)K' \ 

V 2iV, ’ J’ 


s = ,(iV, + l)/2. (15) 


Note that the largest weight oon+i is chosen to be the unpaired one at the center of the 
hfth dimension, such that (5(A) does not change sign for all A E [Xmini^max]- Then the 
approximate sign function S(H) of the effective 4-dimensional Dirac operator (j4]) becomes 


S(H) 


Moll 

1=1 


(h?‘ + C 21 ) 
(K^ + C21-1) ’ 


h 


H 

^min 


(16) 


which is exactly the same as the Zolotarev optimal rational approximation of the sign func¬ 
tion (see Eq. (IT^ ) except for the overall coefficient 


TT (1 + C 2 I- 1 ) 

fj (1 + C2Z) 


(17) 


In Fig. [21 we plot 1 — S'(x) of the optimal DWF with symmetry, for Ng = 17, and 
'Xmax/'Xmin = 6.20/0.05. In Fig. |3](a), the weights are computed according to (IT^ . while in 
Fig. [3](b), according to fiT5|) . In both cases, 1 — «S'(A) has 2?7, -|- 2 = 18 alternate maxima and 
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minima in the interval [Xmin, Amax], with 9 minima and 9 maxima. The maximum deviation 
|1 — S'(A)|max = 2 ^ 2 , where dz is the maximum deviation |1 — y/xR^^’^'^lmax of the Zolotarev 
optimal rational polynomial, as shown in Fig. [H^b). 


IV. NUMERICAL TEST 


Theoretically, the sign function error |1 — S{H) \ of the optimal DWF with symmetry 


is twice of that without R^ symmetry. Nevertheless, it is interesting to chec 
in large-scale simulations of lattice QCD, by computing the residual mass 

1 {tT{[D-\mg)]^{l - S^)D-\mq)}ofl)^ 


c t 


leir difference 


12 | 


TTljrp.pi — 


4r (tr{[(T)c + + mq)-i}o,o) 


u 


(tr(T>c + "^q)o,o> 


u 


m. 


qi 


(tr[75(^c + mq)75(T)c + "ig)]o,o>c/ 

where {Dc + rnq)~^ denotes the valence quark propagator, tr denotes the trace running over 
the color and Dirac indices, and the brackets (•••){/ denote averaging over an ensemble of 
gauge conhgurations. 

To generate the gauge ensemble, we perform the hybrid Monte Carlo simulation of (2-1-1)- 
flavors QCD on the 24^ x 48 lattice with the Wilson gauge action at/3 = 6 / 5 f^ = 6.10 (lattice 
spacing a ~ 0.06 fm), for the sea-quark masses m„a = = 0.005 and rrisa = 0.04, 

with pion mass ~ 260 MeV. For the quark part, we use optimal DWF with R^ symmetry 
for M, d and s quarks, with c = l,d = 0 (i.e., H = Ng = 16, and Xmax/^min = 

6.20/0.05. The weights are computed according to the formula fl9ll. The strange quark is 
simulated with the exact pseudofermion action for one-flavor DWF 10|, and the up and down 
quarks are simulated with the 2-flavors algorithm as outlined in j9|]. We generate the initial 
400 trajectories with 2 Nvidia GTX-TITAN GPUs working together. After discarding the 
initial 240 trajectories for thermalization, we sample one conhguration every 5 trajectories, 
resulting 32 “seed” conhgurations. Then we use these 32 seed conhgurations as the initial 
conhgurations for 32 independent simulations on 32 Nvidia GTX-TITAN GPUs. Each GPU 
generates 100-110 trajectories, and we accumulate a total of 3300 trajectories. From the 
saturation of the binning error of the plaquette, as well as the evolution of the topological 
charge, we estimate the autocorrelation time to be around 10 trajectories. Thus we sample 
one conhguration every 10 trajectories, and obtain 330 conhgurations for this ensemble. 
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Then we compute two sets of quark propagators with point source at the origin and 
rriga = 0.005, one set with the {cus} exactly the same as the sea-quarks (with symmetry), 
and the other set with the {ws} ([2]) without symmetry. The residual mass of the set 
with i ?5 symmetry is {rrireso)^^ = 0.00024(1), while that of the set without R^ symmetry is 
rriresO' = 0.00014(1). Thus the ratio {rrires)^^/'<^res — 1-71(14), consistent with the theoretical 
expectation. 


V. CONCLUDING REMARKS 


With the weights di, dnD and flT^ for optimal DWF with R^ symmetry, together with 
those (E]) for optimal DWF without R 5 symmetry [ll|, this completes the study of DWF 
with the approximate sign function S{H) in the effective 4-dimensional lattice Dirac operator 
satisfying |1 — S'(A)| < 2dz (with R^ symmetry) or |1 — S'(A)| < dz (without R 5 symmetry) 
for A^ G [A^j„, Xmax\y where dz is the maximum deviation |1 — \/x-Rz(3^)1 max of the Zolotarev 
optimal rational polynomial Rz{x) of l/\/x for x G [1, A^^ 3 ,/A^j„], with degrees {n — l,n) 
for Ng = 2 n, and {n, n) for Ng = 2 n + 1 . 

The correspondence between the approximate sign function S{H) of the optimal DWF 
and the Zolotarev optimal rational polynomial gives the optimal rational approximation a 
more general viewpoint. In general, the deviation (5(a:) = 1 — y/xRz{x) can be shifted by a 
constant e, simply by adjusting the overall coefficient Dq in flTSll or do in flT9ll . In particular, 
for if Do is replaced by Dq ([^, then e = +dz, dmax = “^dz and = 0 ; while by Dq 

flTT]) . e = —dz, Smax = 0 and = —2dz- Similarly, for if do is replaced by d'^ flTT)) . 

6 Rdz, dfYiax 2 <dz and Smin O 7 while by d^ fjl4p . 6 d^, d^ax 0 and dYnin 2 dz‘ 

Even though this does not satisfy the criterion that the maxima and the minima of d{x) 
all have the same magnitude and dmin = —dmax, the most salient featues of the optimal 
rational approximation are preserved, namely, the number of alternate maxima and minima 
is (m -|- n -|- 2 ), with {n -|- 1 ) maxima and (m -|- 1 ) minima, as well as all maxima (minma) 
are equal. We can regard this as the generalized optimal rational approximation (with a 
constant shift e). For e > 0, this can be realized by DWF with/without R^ symmetry. 
However, for e < 0 , it cannot be realized by DWF with R^ symmetry and Ng = 2n (even), 
due to the constraint (Ej), (5(A) > 0. 
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Appendix 


In this appendix, we collect the basic formulas of the Zolotarev optimal rational polyno¬ 
mials and for the inverse square root function l/\/a:, x G [!,&]• 


= D„ 


nr.'i'(^+c 


21 


K 




doYl 


nr=i(^+ C21-1 

n. 

x + C 21 

.7 

/=1 


E 

1=1 


B, 


X + C 21-1 ’ 


X + C 21-1 


= [X 


+ C2n) ^ 


1=1 


X + C 2 I-I 


(18) 

(19) 


where 


Ci = 


l-sn2(^;«0’ 


sn 


2! IK' . 


Cl = 


^2n+l 


]k') 


Bj — Dr 


Da = 


u::i{C2.-c2i.i) 

IYi=i,i^iiC 2 i-i - C21-1) ’ 

2^ nr=i(i+^2z-i) 


,2(JJ£_. 

'2n+l’ 
h = dr 


nr=/(c2,-c2z-i) 


nr=i,i7^z(c2i-i - C2z-i)’ 


i+A nr=/(i+c'2z) ’ 


dn — 


2a 


n 


2n 


A=n 




a = 


2n-\-l 

n 


1 + C 21-1 

1 + 1 + C2Z 

(S; ^') 


Here sn(M; k') = rj is the Jacobian elliptic function with modulus k' = — 1/6, as defined 

by the integral 

y/(l-t2)(l-)^’ 

K' = m( 1) is the complete elliptic function of the first kind with modulus k', and 0 is the 
elliptic theta function. 

For the deviation (5(x) = 1 — has 2n -|- 1 alternate change of sign 

in the interval [l,fe], and attains its maxima and minima, 

1 - A 


(5(a:) — dz, —dz, • • • , dz, —dz, dz, dz — 

at consecutive points (1 = Xi < 0:2 < • • • < X 2 n+i = b), 

1 


1 +A’ 


Xi = 


/2.^2 ( (i-m' ■ 


1 — k' sm 


2n ’ 


, i = 1, - ■ ■ ,2n + 1. 


For the deviation 6{x) = 1 — \/xR}^'"'\x) has 2?7, -|- 2 alternate change of sign 

the interval [1,6], and attains its maxima and minima. 


(n,n), 


m 


(5(x) d^, ■ ■ ■ ) dzi dzi dz 


a 


1 + a' 
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at consecutive points {1 = xi < X 2 < ■ ■ ■ < X 2 n +2 = ^), 


Xi = 


1 


1 — k'^su^ 


f {i-l)K' 
2n+l 



1, • • • , ‘In + 2. 
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